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We use the Hopf mapping to construct a magnetic configuration consisting of closed 
field lines, each of which is linked with all the other ones. We obtain in this way a 
solution of the equations of magnetohydrodynamics of an ideal incompressible fluid 
with infinite conductivity, which describes a localized topological soliton. 

1 Introduction 

Solutions of physical equations which have non-trivial topological properties have been stud- 
ied for already more than five years. As examples we may give the "monopole" [1,2] and 
the "instanton" [3] in gauge field theories and the "pseudoparticle" in a two-dimensional 
isotropic ferromagnet [4]. All these solutions are characterized by some topological index: 
the magnetic charge of the monopole and the number of pseudoparticles in the ferromag- 
net are equal to the degree of mapping of a two-dimensional sphere onto a two-dimensional 
sphere, the number of instantons is equal to the Pontryagin index of the mapping of the 
SU(2) group onto the three-dimensional sphere. In each case one can write this index as a 
volume integral of some "topological charge density" . In this connection attention is drawn 
to the integral of motion [5] 



(A is the vector potential) which has been known for a long time in the magnetohydro- 
dynamics of a perfectly conducting fluid and which is called the helicity of the magnetic 
field. Its topological nature is already indicated by the fact that no characteristics of the 
medium in which the magnetic field is present enter into (1). It has also been shown (see 
Refs. [6,7]) that if two field line tubes are linked the integral (1) is proportional to their 
linkage coefficient, i.e., the number of times which one tube is twisted around the other one. 

It is thus clear that the helicity is a topological characteristic of the magnetic field. This 
topological nature of it is completely revealed if we note that (1) is the same as the Whitehead 
integral for the Hopf invariant which characterizes topologically different mappings of the 
three-dimensional sphere 5*^ onto the two-dimensional S*^ [8-10]. The topological meaning 
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(1) 



of the Hopf invariant is simple: it is equal to the linkage coefficient of the curves in 
which are the originals of different points of S^. Hence follows also a more constructive 
conclusion: knowing the mapping — > with a non-zero Hopf invariant, and the simplest 
such mapping was constructed by Hopf himself, we can find the vector field A corresponding 
to it and, then, the magnetic field H = rotA with non-zero helicity. The magnetic field lines 
of this field will be closed and each of them is linked with any other one. In the present 
paper we construct one such magnetic field configuration and we study its properties in 
magnetohydrodynamics. 



2 Stereographic projection and Hopf mapping 

We estabhsh first of all the connection between the physical space and the sphere S^. 
Equation (1) implies that the field H decreases sufficiently fast at infinity so that the helicity 
/ is a gauge-invariant quantity: adding to A the gradient of any function does not change I as 
the additional term after integration by parts gives a surface term which docs not contribute 
because H decreases rapidly, and a volume term which vanishes because divH = 0. If the 
other physical conditions at infinity are also unique (say, we consider a homogeneous isotropic 
medium) we may assume that the Euclidean three-dimensional space is supplemented 
by a point at infinity. Such a "compacted" space becomes topologically equivalent to the 
three-dimensional sphere S^. If we embed in the four-dimensional Euclidean space with 
coordinates (m^; = 1, 2, 3, 4) so that = {u^ : = 1} we can establish the connection 
between R^ and by the stereographic projection 

= — ^, z = 1,2,3. (2) 
I + U4 

It is clear that the point at infinity corresponds to the "south pole" of the sphere with 
coordinates (0, 0, 0, —1). The inverse transformation is reahzed by the formulae 



2xi 1 — X 



^i = T—;^^ "4 = —-^' ^ = 1,2,3, (3) 



1 + x'^' ' 1 + x'^ 

where x'^ — xf is the square of the radius vector. 

Let there now be in R^ a vector field A = {Ai, A2, A^). We find the formulae which 
express the connection between A and the corresponding vector field A^ on S^, where we 
impose on An condition that it be tangent to the sphere: 

u^A^ = 0. (4) 

This condition means that both A and also lie in the tangent spaces to the appropriate 
configuration spaces R^ and S^. We now use the condition for the invariance of the differential 
form ^ 

uja = Adxi = Ai^dUf^. (5) 

^We shall use in this section the convenient formalism of the theory of exterior forms (see, e.g., Refs. [10, 
11]). 
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Since the variables U/^ are connected through the equation = 1 for the sphere in the two 
expressions on the right-hand side of Eq. (5), the number of independent differentials is equal 
to three. Taking as the independent variables on the first three Cartesian coordinates Ui 
we find from (5) 

Ai—^ = Ai + A4—^, ^ = 1,2,3. 



Substituting here (2) and using the equation uf^ — 1, we find three equations 

+ —^^—^A-^A^, 1^1,2,3, (6) 



1 +U4 (1 + U4)^ U4 U4 

which together with (4) are sufficient to express in terms of Ai and vice versa. As a result 

we get ^ ^ 



Ai = {l+U4)A,-UiA4, (8) 

where we must also make the coordinate substitutions (2) and (3) respectively ^. 

We now describe how, using the known mapping f : ^ S"^, one must construct the 
vector field A which corresponds to it and which occurs in Eq. (1) for the Hopf invariant 
of this mapping. One must, as was shown by Whitehead (see Ref. [9] and also Ref. [13]), 
start from the 2-form of the volume on the unit sphere S"^. If that sphere is embedded in 
the three-dimensional Euclidean space with coordinates ^1,^2,^3, the 2-form of the volume 
has the form 

u;2 = (A7r)~\^id^2 A d^3 + ^2d^3 A d^i + ^^^6 A d^) (9) 

(a is the exterior product sign); the coefficient is here chosen in such a way that the integral 
of (9) over the sphere S'^ equals to unity. The mapping / induces a mapping /* in the 
opposite direction from the space of forms on S'^ onto the space of forms on S*^ so that one 
can find the 2-form f*uj2 on S^. One can show that any 2-form on S*'^ can be written in the 
form of an exterior differential of some 1-form uoi where uoi is determined uniquely up to a 
differential of an arbitrary function. We thus can find a form uji such that f*U2 = duji. By 
using a stereographic projection, we can associate with a vector field on S^, determined by 
the form a;i, a vector field A in B? which we can use to evaluate the Hopf invariant through 
Eq. (1). 

We now consider the Hopf mapping f : ^ S"^ which has a Hopf invariant equal to 
unity (see Refs. [8,13]): 

^1 ^2{uiU3 + U2U4), C2^2{u2U3-UiU4), ^3 ^ uj + ul - ul - ul- (10) 

Substituting this formula into (9) we find the form 

f*uj2 — TT~^{dui A du2 — dus A du4), (11) 

^ These formulae can also be obtained from the condition that the covariant derivatives in and are 
the same (Ref. [12]). 
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which is, clearly, the exterior differential of the following form [using the rule for the evalu- 
ation of an exterior differential d{uiduj) = dui A duj = —duj A dui\: 

uji — {2Tr)~^{—U2dui + Uidu2 + 1*4^1*3 — u^du^). (12) 

The vector field corresponding to this form lui = loa [see (5)] 

A = (27r)'"^(-it2, iti, 1*4, -1*3) (13) 

satisfies condition (4) so that we can use Eqs. (8) to find A. Substituting A, thus found, 
into (1), indeed, gives / = 1. However, it is important for us that A, thus found, can be 
identified with the vector potential of a magnetic field with nonzero helicity. 



3 Magnetic field configuration 

We thus find, starting from the vector field (13) and using (8) and (3) a vector potential in 
the three-dimensional physical space. Wc note that Eqs. (3) are clearly written in dimen- 
sionless form, i.e., all coordinates Xi refer to some characteristic dimension R. To change 
to dimensional units we must make the substitution Xi — > xijR^ but in order to keep the 
formulae simple we stay in this section with dimensionless length units. The dimensional 
coefficient of proportionality which fixes the absolute value of the magnetic field strength 
has also so far been dropped. As a result of substituting (13) and (3) into (8) we get, apart 
from a proportionality factor 



_ xxxz - X2 _ X2X3 + xi _ 2x1 + 1- x^ 

'~2(1 + X2)2' '"2(1+^2)2' '~ 4(1 + ^2)2 ■ ^ > 



Calculating the magnetic field corresponding to this potential we find 

_ 2(3:13:3 -X2) _ 2(^23:3 + 3:1) _ 2x^ + 1 - 3:^ 

'~ (1 + 3:2)3 ' '~ (1 + 3:2)3 ' (1 + 3:2)3 ' ^^^^ 

or 

H = rotA/(l + 3:2). (16) 
The square of the magnetic field strength equals to 

SO that the absolute magnitude of the magnetic field of the configuration which we have 
found is spherically symmetric. 

We now find the field lines of the magnetic field (15). The equations of the lines of force 
have the form dx/d/ = H/H, where dl is a line element, or 

dxi _ 2{xiX3 - X2) dx2 _ 2{x2Xs + 3:1) 0^3:3 _ 2x| + 1-3;^ 

"dT" 1 + 3:2 ' "dT" l + 3;2 ' "dT" 1 + 3:2 ' ^ ^ 
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One can easily solve this set of equations if we map it at first on the sphere S^. The vector 
field Hfj^ corresponding to H is found by using Eqs. (7) and (2): 

^1 = -1^2(1 +M4), ^2 = |mi(1 + M4), IU4I + U4), ^2 - -^M3(l + M4), (19) 

SO that the equations for the lines of force on the sphere have the form 
dui du2 dus du4 

d0="^' d0=-"^' ^20) 

where d(j) is the corresponding line element on S^. The solution of the set (20) is clearly: 

= acos(0 + 0i), = sin(0 + 0i), 1^3 = sin(0 + 0o), 1*4 = cos(0 + 0o), (21) 

where the integration constants a and b are connected through the relation + b'^ — 1. 

Again using (2) to change to the physical space we find that the solution of the set (18) 
has the form 

_ acos(0 + 0i) _ asin(0 + 0i) _ 6sin(0 + 0i) 

~ l + 6cos(0 + 0o)' ~ l + 6cos(0 + 0o)' ~ l + 6cos(0 + 0o)' ^ ^ 



where the Z-dependence of (p is found from the differential equation 

d(p 2 



dl 1 + X 



2 



(23) 



which expresses the well-known connection between the line elements in the two metrics: 
the Euclidean and the stereographic (see Refs. [10, 11]). It is clear already from Eqs. (22) 
that the hues of force are closed: when we change from to 27r we completely traverse it 
and return to the initial point. Substituting (22) into (23) we get 

^ = I (24) 

dcj) l + 6cos(0 + 0o)' 

so that the length of a line of force is equal to 



/•2n 

Jo T 



d0 _ 27r _ 27r 
+ bcos(f) ~ (1 - 62)1/2 - y^- 



(25) 



The maximum and minimum values of the radius vector of the points belonging to a line 
of force are found from the formulae 

... I ■ il^'mi'n -1- / /TTf I ... I 



x_-(^^_|^|J , x_-l/x_ ^^^1^1 J . (26) 



The solution of Eq. (24) corresponding to the condition 1(0) — has the form 



a 



a + fa 01 

arctan tan — - — — arctan tan — 

1 + 6 2 / \l + b 2 



(27) 
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whence we find 

tan = ^t^^f + .28) 

a 2 2 

Expressing the trigonometric functions in (22) in terms of tan[(0+0i)/2] and substituting 
(28) we find the way the equations of the line of force depend on /. We shall not write down 
the general formulae in view of their complexity, but restrict ourselves to the case 0i = 0o 
as the curves differing only in the difference 0o ~ 0i can be superposed onto one another by 
a rotation over that angle around the X3-axis: 

(cos 00 + b) cos al — a sin 0o sin al b 
a(l + 6cos0o) a 
(cos 00 + b) sin al + a sin 0o cos al b 

^2 = /, , , > ^3 = -Xl- 

a(i + ocos0oj a 
Hence it is clear that the lines of force are plane curves. Evaluation of their curvature gives 

k=\(fx/df\ = \a\, (30) 

so that the lines of force turn out to be circles of radius l/\a\. This agrees with Eq. (25) for 
their length. 

Although it follows from the way we have constructed the circles that they arc linked, 
it is of interest to verify this also directly. We therefore consider two circles: Ci and C2 
corresponding to values of the parameters a — b — (pi — and different values (po — 
and 00 = 7r/2 (one circle is rotated with respect to the other over 7r/2 around the xa-axis). 
Their parametric equations have the form 



(31) 



The circle Ci lies in the plane X2 = X3, and the circle C2 in the plane Xi = —x^. These 
planes intersect along the line xi — —X2 — —X3. It is clear that if the circles Ci and C2 are 
linked, their points of intersection with this line must alternate with one another. One easily 
finds that Ci intersects this line at the points Ai = (—1, 1, 1) and Bi = (1/3, —1/3, —1/3), 
and C2 in the points A2 = (1/3, —1/3, —1/3) and B2 = (—1, 1, 1). The point A2 lies between 
the points Ai and Bi and the point B2 outside the section {Ai, Bi) so that these pairs of 
points alternate on the line xi — X2 — —X3 and the circles Ci and C2 are linked. 

When the parameter 0o changes from to 7r/2 the circle is shifted in space from Ci to C2 
covering a surface with boundaries Ci and C2 which can be obtained by joining two ends of 
a strip after twisting it over 360°. It is known (and one can easily verify this experimentally) 
that if one cuts such a strip along its boundaries following a closed line it falls apart into two 
such strips which are linked. Continuing this cutting exercise we shall obtain ever narrower 
strips which are linked with one another. It thus becomes clear that all circles forming the 
original strip arc linked with one another. 

When the parameter 0o changes from to 27r the circle describes a closed surface (a torus 
obtained from a cylinder which is twisted 360° before it ends are joined) which is bounding 



Ci = jxi = \/2 cos(//a/2) — 1, X2 = sin(//A/2), 0:3 = sin 
C2 = |xi = -sin(Z/V2), X2 = V2cos(VV^) - 1, X3 = sm{l/V2 
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a "plait" of closed lines of force. The lines of force thus lie on toroidal surfaces which are 
imbedded one into another, and are circles, each of which is linked with all the others. 

We now consider a physical system in which the magnetic field configuration which we 
have described can be realized. 



4 Magnet ohydro dynamic soliton 

We change to dimensional units so that the magnetic field (15) takes the form 

H = (^^^{2^[k X x] + 2(k • x)x + {R' - x')k}, (32) 

where k is the unit vector along the a;3-axis, R the size of the sohton, and Hq the magnetic 
field strength at the origin. The square of the magnetic field strength is equal to 

H^R^ 

^ (R-lx-r ^^^^ 

Using Eq. (1) to evaluate the helicity of the magnetic field we get 

I = ^H',R\ (34) 

We note that through the mapping x — x we get an "antisoliton" , the magnetic field of 
which differs from (32) in the sign in front of the first term in the braces, while the helicity 
(34) also changes sign. 

We shall consider a perfectly conducting liquid for which I is an integral of motion. 
We also restrict ourselves to the case of an incompressible ideal fluid. The equations of 
magnetohydrodynamics for stationary flow have the form (see, e.g., Ref. [14]) 

divH = 0, divv = 0, rot [v x H] = 0, 

They are clearly satisfied (see Ref. [15]) when the fluid moves along the magnetic field lines 
of force with a velocity 

while the pressure satisfies the equation 

P + — = Poo = const. (37) 

OTT 

Thus, Eqs. (32), (36), and (37) give an exact solution of the equations of magnetohydrody- 
namics which describes a localized topological soliton. 
We evaluate the soliton energy 



7 



For a physical interpretation of topological solitons we must bear in mind that they are 
metastable states, the energy of which is higher than the energy of a state at complete 
equilibrium. It is thus necessary for the stability of a soliton, at any rate, that there does 
not exist such a continuous deformation at which its energy diminishes while the topological 
invariant is conserved. Comparison of (38) and (34) shows that 

so that the soliton can diminish its energy for constant / by increasing its radius. However, 
in the case considered there is yet another integral of motion — ^the angular momentum [we 
take the + sign in Eq. (36)] 



M = p y"[x X w]d^x = ^{p7r^y/^HoR% 



which stabilizes the "spreading" of the soliton (cf. the remarks about "collapse" of solitons 
in Refs. [16, 17]). 

The radius R and the field Hq are completely determined by the two conserved quantities 
I and M: 

R^{My4npiy/\ Ho^8{p/TTy/^I/M. (39) 

For given /, the M-dependence of the energy has a specific decreasing spectrum, E oc M~^/^. 
One must, however, bear in mind that / and M are not completely independent quantities. 
As the pressure is always positive, it follows from (37) and (33) that 

< 87rpoc (40) 
and thus, according to (39), / and M must satisfy the thermodynamic inequality 

T / \ 1/2 

/ TT / Poo ^ 



— < - — • (41) 

For a given external pressure Poo the radius and energy of the soliton satisfy thus the 
inequalities 

R > (27^Voo)-'/'M^/^ E < {nyj2''py/'M'/\ (42) 
Combining these inequalities [or substituting (40) into (38)] gives 

E < {7i/2fp^R^ (43) 

which is essentially the same as the well known inequality E < 3pV which follows from the 
fact that the trace of the energy-momentum tensor is positive (see Refs. [18,19]). 

The magnetic field of the soliton (32) is produced by currents which circulate along closed 
lines with a density 

j = ^rotH = — — ^— (2i?H + H X x). 

Att 2'kR^ + x^^ ' 

These currents are conserved since we neglect dissipative processes. When account is taken 
of the finite conductivity cr, magnetic field diffusion occurs. The considerations given here 
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are applicable if the hydrodynamic velocities dominate the diffusion velocities, i.e., when 
{ujnR/vy/'^ <^ R, Ujn — c^/4:7ra is the magnetic viscosity, or 



Re, 



■m 




(44) 



the magnetic Reynolds number must be much larger than unity. When this criterion is 
satisfied, the condition aR/v 3> 1 that the displacement current is negligible (sec Rcf. [14]), 
which is assumed to be true in magnctohydrodynamics, is satisfied automatically (the dis- 
placement current vanishes identically in a stationary case when there is no dissipation). We 
can estimate the lifetime of the soliton by dividing its energy E by 



As a result we get 

t ~ — . 45 

When applying inequality (44) to this problem this means that the lifetime (45) is much 
longer than the characteristic time ^ R/v for the motion of a fluid particle along a line of 
force. 

5 Conclusion 

The equations of magnetohydrodynamics thus admit of an exact solution which describes a 
localized topological soliton. This kind of solution has already been met with in the physics 
of the condensed state (see, e.g., Refs. [16,17]). We note here some difference between the 
magnetohydrodynamic soliton and, say, a soliton in a ferromagnet [17]. In a ferromagnct the 
mapping — > S"^ is realized by the order parameter — the magnetization vector m.{x). Here 
the sphere S"^ has a direct physical meaning, namely, it is the configuration space of the vector 
m. At the same time the map of a point from S"^ has no special physical meaning — it is the 
line on which m(x) takes a constant value and the Hopf invariant characterizes the linking 
of such lines. In magnetohydrodynamics there is no ordering parameter and the sphere S'^ 
has a completely arbitrary character: its points merely "number" the magnetic lines of force 
and the correspondence between the hues of force and the points on S'^ is established by the 
Hopf mapping S"^ . This mapping is not realized in such an apparent manner as in 

the case of a soliton in a ferromagnet, but now the maps of the points of 5*^ have a direct 
physical meaning — they are the magnetic lines of force, and the Hopf invariant characterizes 
their linking. 

I express my gratitude to A. A. Vedenov, V.G Nosov, A.L. Chernyakov, V.R. Chechetkin, 
and V.V. Yan'kov for discussions of the results of this paper. 
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